Phase-transitions in isotropic extreme type-II superconductors 
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Using large scale Monte Carlo simulations on a uniformly frustrated 3DXY model, we report a first 
order vortex lattice melting transition in clean, isotropic extreme type-II (ft —> oo) superconductors. 
This work clarifies an important issue: the unpinned vortex liquid is always incoherent with no phase 
coherence in any directions for all anisotropics. Previous claims of a disentangled vortex liquid for 
isotropic superconductors based on simulations, are due to finite size effects. We explicitly show that 
the effective vortex-line tension vanishes precisely at the superconducting phase transition in zero 
magnetic field. This loss of line tension is accompanied by an abrupt change in the connectivity of 
the vortex tangle across the superconductor. We also obtain results indicating that the connectivity 
of the vortex tangle changes in a similar way even in finite magnetic field, and suggest that this 
could also be associated with a genuine phase-transition. 



I. INTRODUCTION 

Since the discovery of copper oxide based high tem- 
perature superconductors (HTSC)El, which are of the 
extreme type-II variety, there has been great interest 
in their phenomenolo&ical phase diagrams. Abrikosov's 
mean-field descriptiono, which is valid for conventional 
low temperature superconductors, is expected to be mod- 
ified by the strong thermal fluctuations in HTSC. Exten- 
sive research in both theory, numerical simulations, and 
experiments over the years have resulted in a general con- 
sensus on some of the fundamental issues. The current 
understanding of HTSC in a uniform magnetic field B is 
as follows. In the absence of any pinning disorder, the low 
temperature Abrikosov vortex lattice phase melts into a 
vortex liquid via a first order transition at the tempera- 
ture T m . The Abrikosov vortex lattice phase is charac- 
terized by a transverse triangular crystalline order and a 
finite longitudinal phase coherence. 

However, recently there has been some debate about 
the nature of the vortex liquid which the Abrikosov vor- 
tex lattice melts into as temperature iacxeases. Numer- 
ous simulations using the 3DXY modetTQ, LondonQ and 
Lowest Landau Level approximationsu have indicated 
that the vortex liquid is incoherent, i.e. the phase coher- 
ence or superfluid density in any direction, is zero. The 
crystalline order and phase coherence are destroyed si- 
multaneously at the melting transition. This scenario has 
been supported by experiments on high quality YBCO, 
crystalsa. Other simulations using the 3DXY modelEHliLHl 
have suggested that the longitudinal phase coherence per- 
sists above the melting transition and only vanishes at a 
higher "entanglement" temperature Te- In this scenario, 
the vortex liquid at T m < T < Te would be disentan- 
gled with relatively straight vortex lines. For T > Te, it 
was claimed that larger thermal fluctuations causes the 
vortex lines to be entangled with a concomitant loss of 
global phase coherencet3~li3. 

Previously, a vortex liquid with nonzero longitudinal 



phase-coherence, or superfluid density, was found irusim- 
ulations on an isotropic systemE£rE3. More recentlyE3, it 
has been proposed that for large B, T m and Te merge 
into a single transition whereas for small B they are well 
separated. The authors found that for an isotropic sys- 
tem, the crossover between the two regimes is at a filling 
fraction 1/18 < / < 1/6, where the precise definition of 
/ oc B will be given shortly. However, these simulations 
were performed on rather small systems. In the present 
work, we have performed similar simulations on a much 
larger system, and found an exclusively incoherent vortex 
liquid down to / = 1/60. This implies that the only ther- 
modynamically stable vortex liquid phase, is one which 
has zero longitudinal superfluid density, with full trans- 
lational and rotational symmetries and zero phase coher- 
ence in all directions. We are therefore led to the conclu- 
sion that the recently discussedtj crossover between the 
phase-coherent vortex liquid and the phase-incoherent 
vortex liquid is a numerical artifact. 

The organization of this paper is as follows. In Section 
II we first briefly review the model used in the simula- 
tions, and the quantities we calculate. In Section III wc 
discuss our results and their implications. In Section IV 
we provide a summary and conclusion. 



II. MODEL 

We use a uniformly frustrated 3DXY model on a cubic 
lattice to describe an isotropic, extreme type-II supercon- 
ductors in a magnetic field. The London model for super- 
conductors can be readily derived from the phenomeno- 
logical Ginzburg-Landau (GL) model with the approx- 
imation that the amplitude of the local complex order 
parameter ip = \tp\exp[iff\, is fixed. The Hamiltonian H 
of the London model consists of degrees of freedom in the 
phase 6(r) and the gauge vector potential A„ p (r) asso- 
ciated with the magnetic induction B inside the system, 
i.e. H = H[{6(r), A vp (r)}}. 



1 



For an isotropic extreme type-II superconductor, the 
penetration depth A is much larger than the coherence 
length £ such that the GL parameter k — A/£ — > oo. This 
means that the magnetic fields surrounding the vortex 
lines strongly overlap with one another giving a spatially 
smooth B. This condition is ensured in the regime where 
B > B c i. In other words, the fluctuations of A„ p (r) on 
the length scale of £ are negligible compared to the fluc- 
tuations of 8(r). Therefore, we can further simplify the 
problem by dropping the degrees of freedom of A vp (r) 
from the Hamiltonian H, and fix B equal to the exter- 
nal applied magnetic field. The resulting Hamiltonian 
H = H[{6(r)}} is the 3DXY model. This k -► oo (or 
frozen gauge) approximation has been widely used as a 
phenomenological model for superconductors. Note that, 
within this model, the system has no magnetic flux lines, 
since there can be no tubes of confined magnetic flux 
when A — > oo. The system only exhibits vortex lines. 

In order to perform numerical simulations on the re- 
sulting model, we discretize the model on a 3D cubic lat- 
tice with grid spacing £. The di|B^<jnsi,onlcss Hamiltonian 
of this lattice model is given by 



H[{6(r)}] 



-Jo 



E 



(1) 



where Jo = ( f>Q^/47r 2 /ioA 2 is the isotropic coupling energy, 
<&o = h/K / e is the flux quanta and fj,o is the permeabil- 
ity of a vacuum in SI units. The dimensionless vector r 
labels the position of an arbitrary grid point. We define 
the gauge-invariant phase difference j a (r) as 



j a (r) = 9(r + e a ) - 9(r) - 



2tt£ 



r+e Q 



rfr'-A^(r') (2) 



where e Q is the unit vector along the a-axis. The con- 
vention is that ja(r) "flows" from the grid at r to r + e Q . 
In this model, it is natural to define a dimensionless tem- 
perature T = ksT/J{). The magnetic induction B is 
conveniently represented by the filling fraction 
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In this paper, we perform a Monte Carlo simulation 
on an isotropic system with V — N x x N y x N z num- 
ber of grid points. The aim of our numerical simulations 
is to identify possible thermodynamic phases (i.e. vor- 
tex lattice or liquid, etc.) and phase transitions (first 
order or continuous) or crossovers associated with them. 
Specifically, we are interested in calculating the internal 
energy, specific heat, structure factor, and helicity mod- 
uli. These thermodynamic quantities and their physical 
significance will be discussed next. In addition, we in- 
troduce and define a quantity Ol which we denote the 
vortex-path probability, and discuss some of its physical 
implications. 

A first order transition is indicated by a (5-function 
anomaly in the specific heat, equivalently a discontinu- 
ous jump in the internal energy. On the other hand, the 



hallmark of a continuous transition is a jump in the spe- 
cific heat, modified by fluctuation contributions to the 
anomaly, and a continuous internal energy. The inter- 
nal energy per site is obtained by averaging H in Eq. (Q) 
over the thermal equilibrium states, normalized by the 
total number of grids, i.e. E =< H > /V. We define 
a dimensionless specific heat per site C using the stan- 
dard fluctuation theorem of a classical system with Gibbs 
distributionEJ: 



C = 



< H 2 > - < H 

v{k B Ty 



(4) 



A convenient and widely used quantity to probe the 
global phase coherence of the system is the helicity mod- 
ulus, which is proportional to the second derivative of 
the free energy associated ptiih Eq. (|f|) with respect to 
an infinitesimal phase twisttd'S. On a square lattice, the 
dimensionless helicity modulus T Q can be written as 
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Tqi measures the "stiffness" of 9 with respect to twist- 
ing along the a-direction. If T a > then 9 is stiff in the 
a-direction, or more precisely, there is global phase coher- 
ence or superconducting response in that direction. By 
tracking the temperature at which T goes to zero, one can 
determine the superconducting-normal metal transition 
temperature of the system. In the mixed phase T|| > 
and = 0, where the subscripts || and _L denote the 
directions longitudinal and transverse to B, respectively. 

The structure factor probes the transverse crystalline 
order of the vortex system. We adopt a conventional 
definitions: 



S(k ± ) 



1 
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^ v\\ (r) exp(ikj_ • r) 
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where kj_ is a two-dimensional reciprocal vector and 
v a (r) is the,k>cal vorticity measured on the dual square 
lattice gricoo, composed of the centers of every direct 
unit cell. A crystalline ordered phase is characterized by 
S'(ko) > (or Bragg peaks) where ko are the discrete set 
of reciprocal lattice vectors associated with the crystal 
structure. On the other hand, 5(k) for a phase with full 
rotational invariance exhibits ring patterns. 

At a fixed temperature, the equilibrium configurations 
are generated by making random changes to (9(r) at each 
grid-point via the the Metropolis algorithm. This is 
equivalent to randomly changing all the six j a attached 
to each grid-point. To ensure conservation of vorticity 
in each unit cell, j a has to be in the range [—it, 7r). Ad- 
dition of ±2-7r shall be administrated to bring j a back 
into range at every Monte Carlo step. This procedure 
introduces vortex-loops into the system. Such loops are 
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the elementary topological excitation of the model. Peri- 
odic boundary conditions (PBC) are imposed on j a such 
that j a (r + N^e^) — j a ( r ) f° r M = (^i V, z )- More details 
about the Monte Carlo procedure may be found in Refs. 

IS 

We define Oi as the probability of finding a directed 
vortex path threading the entire system transverse to 
the induction B, without using the PBC along the field 
directiorilB. It is obtained by computing the number 
Ny of times we find at least one such path threading 
the system in any direction _L B in Np different phase- 
configurations, normalized by Np, i.e. Ol = Ny/Np. 
The fact that Ol = implies that there is no con- 
nected path of vortex segments that threads the entire 
system in the transverse direction, without using PBC 
along the field direction several times. Now, let N" 
(a G [x,y, z\) denote the areal density of connected vor- 
tex paths threading the system in any direction, includ- 
ing the direction parallel to the induction. It is clear 
that in the Abrikosov vortex lattice phase Ol = 0, and 
N[ = B/$ , while Nl = N V L = 0. Thus, N£ is a con- 
served quantity at fixed induction B. On the other hand, 
Ol = 1 implies that N^ ,y > 0, and the total number 
of vortex paths threading the system in any direction 
scales with system size, but undergoes thermal fluctua- 
tions. Therefore, N" is no longer a conserved quantity. 

We propose the following scenario to interpret 
the change in Ol and N". Number conservation 
uniquely identifies a [/(l)-symmetry, and hence the low- 
temperature phase of the vortex-system (the dual of the 
phase-representation of the superconductor) exhibits ex- 
plicit [/(l)-symmetry, since Ol — 0. At high temper- 
atures Ol = 1, is not conserved, and the U(l)- 
symmetry is broken. A [/(l)-symmetric phase cannot 
be analytically continued to a [/(l)-nonsymmctric one. 
The change in Ol from to 1 could therefore signal a 
phase-transition, in this case involving breaking a global 
{/(l)-symmetry, in finite as well as zero magnetic field. 
However, to substantiate such a claim, one needs to ar- 
gue that Ol is related to a local order parameter of the 
system. To this end, we note that it is possible to tran- 
scribe the vortex-part of the Ginzburg-Landay theory in 
the phase-only approximation in such a way that the 
vortex-part of the theory is specified in terms of a local, 
complex matter field </>(x), and that the theory then ex- 
plicitly exhibits a L/(l)-symmetryE3. In the Lattice Lon- 
don Model, corresponding to the Villain-approximation 
to the Ginzburg-Landau theory, this symmetry is there- 
fore only implicit, or "hidden" . The probability of find- 
ing a connected vortex-path starting at a point x and 
ending at point y, G(x, y), is given in terms of the 
two-point correlation function of the matter field ^>(x), 
G(x,y) =< </>*(x) (j)(y) >t£l. Ol may be viewed as 
a special case of G, and in the thermodynamic limit 
corresponds to limix-yi^oo G(x, y). If G(x, y) is non- 
zero in this limit, this suggests the possibility of hav- 
ing < <f>{x.) >^ 0, and hence a broken t/(l)-symnietry. 
Although this does not constitute a proof that Ol is 



connected to a local order parameter whose expecta- 
tion value is associated with a broken t/(l)-symmetry, it 
seems to be suggestive of such a phase-transition existing 
even at finite magnetic field. Note that the above local 
matter field cj> appears to be the dual field of a complex 
order parameter appearing in a somewhat different inde- 
pendent approach to the same problem by TesanoviuLZL 

III. RESULTS AND DISCUSSION 

In this section, we discuss the results of our simulations 
on an isotropic system with size of 120 3 grid points, which 
is the largest to date on an isotropic system. The system 
is subdivided into multiple sections, and the Monte Carlo 
procedure is implemented in parallel across, the sections 
using 3D "black and white" decompositiontl The filling 
fraction considered are 1/f = 20,40,60, and oo. The 
system is cooled from high temperatures. For each tem- 
perature, a typical run consists of 120000 Monte Carlo 
sweeps across the whole lattice, 1/3 of that were used 
for equilibration. Near phase-transitions, up to 600000 
sweeps were used. 

In Fig. (|l|), we present the results for the zero-field 
case, / = 0, where we show the specific heat G, the 
superfluid stiffness T, and the vortex-path probability 
Ol as functions of temperature. Note how the specific 
heat anomaly, the vanishing of the superfluid stiffness 
and Ol, all coincide with T c 2.20. The physical inter- 
pretation of Ol = 1 is that Ihxj-cffcctive long- wavelength 
vortex-line tension vanishesa&EZl. This claim is substan- 
tiated by calculating the probability distribution D(p) 
of vortex loops as a function of perimeter p, at various 
temperatures. We may fit this distribution to the formES 

D{p) = A p- a e- pep ; T < T c 

= A p- a ; T>T C (7) 

where s is an effective, temperature dependent, long- 
wavelength line tension, p = 1/kpT, and a — 5/2. 
The results are shown in Fig. (|2|), demonstrating that at 
T = T c , a purely algebraic decay is realized, implying 
e = 0. Below T c , the exponential decay is well fitted, 
showing that e ^ 0. The inset of Fig. (|^) shows e as a 
function of temperature, as obtained in our simulations. 
A similar method of extracting e for the zero-field case 
using a similar form for D(p) (without he power-law pref- 
actor) has previously been used by Li and Teitel, in Ref. 

The results for the quantities S(k), T and G for 
1/f = 20 are shown in Fig. (^). Similar qualitative fea- 
tures are also found for the cases 1/f — 40,60. First, 
<5(kj_) exhibits six-fold Bragg peaks at low temperatures 
(not shown). This indicates that the low temperature 
phase is a triangular vortex lattice phase. The destruc- 
tion of the vortex lattice structure is marked by a melt- 
ing transition at T m rj 1.34 where <S'(ki) drops sharply 
to zero, where ki is the wavevector of one of the six first 
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order Bragg peaks. The sharpness of the drop in S'(ki) 
strongly suggests that the transition is first order. This is 
confirmed by the appearance of a J-function like peak in 
C at the same temperature T m . Coincidentally, T|| which 
is finite at low temperatures, also drops sharply to zero at 
T m . The isotropic system exhibits longitudinal supercon- 
ductivity below T m , but not above. Moreover, based on 
our above discussion in Section II, we interpret the rise in 
the quantity Ol from to 1 as signalling that the effec- 
tive long-wavelength vortex-line tension vanishes. Thus, 
the vortex liquid phase is divided into two regions in 
phase-space. In one region, the vortex liquid is phase- 
incoherent i.e. has no longitudinal superfluid density but 
has finite vortex-line tension. We propose that this phase 
exhibits a vortex-associated t/(l)-symmetry. In the other 
region the vortex liquid is phase-incoherent but has zero 
vortex-line tension. We propose that this phase exhibits 
a broken [/(l)-symmetry. 

Based on this, one can conclude that a first order melt- 
ing transition of the Abrikosov vortex lattice exists in 
an isotropic system in the absence of pinning. The en- 
tire vortex liquid phase is incoherent, i.e. a vortex liquid 
phase with no longitudinal superconductivity. Note that 
the same conclusion has been. reached in earlier simula- 
tions on anisotropic systemsouu. There have been earlier 
reports of a disentangled vortex liquid, i.e. vortex liquid 
with non-zero longitudinal phase coherence. jEfrfse re- 
sults have been obtained in similar simulationll3il3. In 
Rcf . [l^, it was argued that a phase-coherent vortex liquid 
should be most easily observed in isotropic systems for 
/ < 1/18. However, these simulations were performed on 
comparatively small systems, typically of sizes no larger 
than 24 3 . Our present results are based on much larger 
systems, 120 3 . We have observed a slight difference in 
the temperatures at which S(k) and T|| vanish in sim- 
ulations on smaller system or lower number of sweeps. 
Therefore, it may be concluded that the existence of a 
vortex liquid with a nonzero longitudinal superfluid den- 
sity is a numerical artifact _af small system sizes and/or 
insufficient simulation timeEj. It is well understood that 
the A-transition is driven by proliferation of thermally 
excited loops of all sizes. In recent work by some of us, 
it was proposed that the same mechanism is driving the 
first order melting transition at low B and a newly dis- 
covered continuous transition involving the breaking of 
U(l) symmetry at large Bu. 

Similar results for l/f = 40 and 60 enable us to pro- 
pose a simple phase diagram for an isotropic extreme 
type-II superconductor in the absence of disorder (see 
Fig. (||)). Immediately below and above the T m line, 
the phases of the vortex system are identified as the 
Abrikosov vortex lattice and the incoherent vortex liq- 
uid, respectively. 

In a pin-free system, one would expect Tj_ to be zero 
at all temperature. In this case, the numerical lattice, on 
which simulations are performed, effectively pins the vor- 
tex lines from moving in the transverse plane and coun- 



teracts the Lorentz forces on them in the presence of a 
transverse applied current. However, the pinning is over- 
come by thermal fluctuations at higher temperature and 
the vortex lines are depinned at temperature Td w 0.62. 
Fortunately, we see that Td <C T m , which means that 
near T m , the vortex lines, and the melting process, are 
completely free from the effects of the numerical grid. 
Therefore, the features of C, T|| and S(k) at T m are gen- 
uine thermodynamic effects. 



IV. CONCLUSION 

We have performed simulations of the uniformly frus- 
trated 3DXY model on a large isotropic system (120 3 grid 
points) for a variety of filling fractions l/f — 20,40,60 
and oo. We found a first order melting transition in this 
isotropic system for all the three non-zero values of / con- 
sidered. The longitudinal phase coherence and triangular 
crystalline order of the Abrikosov vortex crystal are si- 
multaneously destroyed at the melting transition. Above 
the melting temperature, the incoherent vortex liquid is 
the only thermodynamic phase. We have demonstrated 
that previous-claims of the existence of disentangled vor- 
tex Iiquidll2nl3 is due to performing simulations using in- 
sufficient system sizes and simulation times. 

We have shown that the effective vortex-line tension 
vanishes precisely at the zero-field superconducting tran- 
sition. The loss of superfluid stiffness, the loss of line 
tension, and the abrupt change in the connectivity of the 
vortex tangle, as signalled by the change in the quantitity 
Ol across the system, all coincide in this case. A simi- 
lar change in connectivity across the vortex system takes 
place at finite magnetic field. The results of the Ref. 
|l5| , the present paper, and in particular those of Ref. ^ 
strongly indicate that this change in connectivity is sharp 
in the limit of large systems, thus indicating the loss of 
number-conservation of connected vortex paths thread- 
ing the system. Since the finite probability of finding a 
connected vortex path threading the system in a direction 
other than the magnetic field may be tied to the-.finite 
expectation value of a local complex matter fielcfij, this 
lends further support to the argument that the change 
in Ol signals the breaking of a [/(l)-symmctryEJ. At 
the very least this proposition appears to be intriguing 
enough to warrant further investigation. 

We finally caution the reader that we so far have not 
been able to detect any anomaly in specific heat at the 
suggested new finite-field transition inside the vortex liq- 
uid, for the isotropic case. Even the anomaly at the first 
order melting transition is weak in the isotropic case, and 
is difficult to bring out in simulations. It may be that 
considerably larger systems are needed for the isotropic 
case in order to see signals in the specific heat of the 
suggested new transition due to the small amount of 
entropy in the transition. This is the reason why the 
anomaly at the first order vortex lattice melting transi- 
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tion is difficult to observe in simulations. It is conceiv- 
able that increasing the anisotropy of the system should 
bring out the anomaly clearer, if it exists. This indeed 
is the case for the anomaly at the vortex lattice melting 
transition. A weak anomaly associated with the putative 
?7(l)-transition, may in fact have been observed for the 
anisotropic case in Ref. [| 
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FIG. 1. Specific heat C, superfluid stiffness T, and vortex-path probability Ol for T = 1, / = 0. 



FIG. 2. Vortex-loop probability distribution D(p) as a function of loop-perimeter p for various temperatures, F = 1, / = 0. 
The lines in the figure are fits using -D(p) = A p~ a exp(—f3ep), with A — 1, a — 5/2, and /3 = 1/fceT. e is the only fitting 
parameter in all plots. The inset shows the effective long-wavelength vortex-line tension e. The solid line is a guide to the eye. 
e vanishes at T = T c . Here, £ is the grid-spacing of the numerical lattice, and serves both as a unit of length and a measure of 
the superconducting coherence length. 

FIG. 3. Specific heat per site C, helicity moduli Ti and |, vortex-path probability Ol, and structure factor 5(ki) (where 
ki is the wavevector for one of the first order Bragg peak) as a function of T for a system size V = 120 3 and 1// = 20. The 
melting temperature T m « 1.34 is marked by the sharp drop of S'(ki). The coincidence of a sharp peak in C at T m confirms 
that the melting phase transition is first order. Tii also vanishes at T m indicating that the triangular vortex crystal melts 
into an incoherent vortex liquid. At T Ri 1.90, Ol jumps from to 1, signalling a [/(l)-symmetry breaking, or equivalently a 
3-DXy-transition. 



FIG. 4. The intrinsic f — T phase diagram of an isotropic system based on simulations on system size V = 120 3 and 
/ = 0, 1/20, 1/40 and 1/60. The first order melting-line and the 3DXY-line are denoted by T m and Ti, respectively. 
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